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Abstract. This note deals with concrete applications of the general treatment of anisotropic 
thcrmo-elasticity developed in the first part [M. Reissig, J. Wirth, Anisotropic thermo- elasticity 
in 2D - A unified treatment, Asympt. Anal. ?? (????) ??—??]. We give dispersive decay rates 
for solutions to the type-1 system of thermo-elasticity for certain types of anisotropic media. 

r- 
o 
o 

J> 1. Elastic operators and thermo-elastic systems 

o 

In the first part, [5] we treated the type-1 system of thermo-elasticity for an anisotropic medium 
in two space dimensions. This system reads as 

^ u tt + A(D)U + jV6 = (1.1a) 

6 t - kA6 + 7 V • U t = (l.lb) 

for the elastic displacement U(t, x) € R 2 and the temperature difference to the equilibrium state 
9(t,x) € R. Physical properties of the medium are described by the thermal conductivity k > 0, 
i-Q the thermo-elastic coupling 7 > and the elastic operator A(D). Here we assume that its symbol 

& has the structure 

A(£)=D(0 T SD(£), (1.2) 

where 

CM Ai , 

> D= 6 I (1-3) 

is a matrix of first order symbols and 

fn A 

OO S = A T 2 (T 2 I (1.4) 

\CTl cr 2 

contains the elasticity modules of the medium. Usually one makes the assumption that S is 
positive definite such that A(D) is a positive and self-adjoint operator. Then the first equation 
is hyperbolic, while the second one is parabolic. 

In [2] the thermo-elastic system was treated under the assumptions 

(Al-2) A(£) is positive (self-adjoint, 2-homogeneous and real-analytic in 77 = £/|£|); 

(A3) #spec A(t)) = 2 for all 77 e S 1 ; 

(A4) 7 2 7^ Xj(v) — tr -4.(77) f° r a H hyperbolic directions (with respect to ^(77) € spec A(rj)), 

where a direction 77 e S 1 is called hyperbolic with respect to Xj(rj) if the corresponding (nor- 
malised) eigenvector 7^(77) is perpendicular to the direction 77. Under these conditions the system 
was reformulated as system of first order and Fourier integral representations in terms of Fourier 
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integral operators with complex phases established. They in turn lead to L p ~L q decay estimates 
of solutions. 

Here we will focus on these results. It turned out that the decay properties of solutions depend 
only on the neighbourhood of hyperbolic directions. Essential are the vanishing orders tj(rj) 
of the so-called coupling functions a.j{rj) — rj ■ Tj{rj) for the corresponding eigenvectors and the 
orders of tangency jj (77) for the Fresnel curves 

S = { £ | 1 € spec A($ } = { ^(r?)" 1 / 2 ?? | rj G S 1 , xjfa) € spec A( V ) } (1.5) 

in the hyperbolic direction and for the corresponding sheet. The following theorem collects 
dispersive decay rates, the statements in [5] are more localised in the sense that we decomposed 
solutions into modes with different decay behaviour micro-locally. 

Theorem (Reissig-Wirth, [5]). Let fj £ S 1 be a fixed direction, x G Co°(^ 1 ) a cutoff function 
for a small neighbourhood of fj and = x(£/l£D- Then the following micro-local dispersive 
estimates hold true: 

(1) If fj is hyperbolic with respect to Kj{fj) and 2£j(fj) < Jj(fj), then for any r > 2 we find a 
constant C such that 



|| X (D)([/ f , y/A^)U, 0)(t, Olloo < C(l + t) 2 ^> (||x(D)(C/ f , ^AlpjU, 9)(0, -)\\ m ,r). (1.6) 

(2) If fj is hyperbolic with respect to Kj{fj) and 2tj{rj) > Jj(fj), then for any r > 2 we find a 
constant C such that 



\\xp)(U t ,y/A{piU,9)(t, ■)]]<*> <C(l + t) ^ (\\x(T>)(UuVA^)U,d)(0,-)\\ H i, r ). (1.7) 
(3) If fj is not hyperbolic, then for any r > 2 we find a constant C such that 

||x(D)(£/ t , y/A(pjU, 9)(t, OHoo < C(l + t) -1 (||x(D)(f/t, y/A(p)U, 0)(O, (1.8) 

We want to apply this statement in several concrete situations. We prepare this in Section [2] 
by some general remarks related to the choice of (1.2) and the resulting achievable decay rates. 
Later in Sections [3] to [6] we will discuss isotropic, cubic, rhombic and a case of fully anisotropic 
media. 



2. Some general statements 

The Fresnel curve S has the form S — Si U S 2 with Sj = {xj(ry) _1 / 2 77 | rj e S 1 } and due to (A3) 
we are allowed to label the (real-analytic) eigenvalues in ascending order, H\{r\) < Hi{r\). 

Because the entries of the matrix A(£) are polynomial of degree two we immediately see that the 
associated Fresnel curve S is algebraic of degree four. Hence, any straight line intersects S in at 
most four points. As pointed out in pP this has important consequences. Any line intersecting 
the inner sheet 5*2 has to intersect the outer one S\ in two points. Thus it can intersect S2 at 
most twice and therefore we see that S2 has to be strictly convex. Thus the order of tangency 
72(77) for tangents on S2 at ^(rj) -1 / 2 ^ is two for any direction 77. Similarly we see that for the 
outer one 2 < 71(77) < 4. 

The hyperbolic directions are determined by the fact that r\ and i] 1 - are eigenvectors of the 
matrix A(rf), thus they satisfy rj 1 - ■ A{jj)rj = 0. This is a polynomial equation of order eight. 
Introducing polar co-ordinates 77 = (cos </>, sin </>) T allows to rewrite it as trigonometric polynomial 
and hyperbolic directions are zeros of 

4(<7i + <t 2 ) cos 20 + 4(<7i - <j 2 ) cos 40 - 2(ti - t 2 ) sin 20 - (n + r 2 - 2A - 4/x) sin 40. (2.1) 
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It turns out that this polynomial has either four, six or eight roots according to the choice of S 
or it vanishes identically (for isotropic media). 



The vanishing orders of the coupling functions are bounded by the vanishing order of (2.1|: If 
a>j{rj) vanishes to order k, then Tj{rj) = rj up to order k and therefore 77 • A^rj 1 - vanishes at 



least to order k. But this means that the order of zeros to (2.1 ) always dominates the vanishing 
order of the coupling functions and therefore (provided that the medium is anisotropic) 

4 ^ E E *M<* ( 2 - 2 ) 

j=i,2 77:03 (n)=o 
for the vanishing orders £j {rj) of ctj in rj. 

These general observations imply that Theorem [T] gives as overall decay rate for the medium 
either i~ 1//2 , t -1 / 3 or £ -1 / 4 . We will provide examples for all these decay rates. 



3. Isotropic media 
For isotropic media we have a\ = 02 = and t\ = t 2 = A + 2/i, i.e. 

In this case spec A{rj) = {A + 2ji, ji} with corresponding eigenvectors rj and 77- 1 -. Thus one of the 
coupling functions (the one corresponding to /1) vanishes identically and the Fresnel curves are 
two concentric circles with tangency indices 71(77) = 72(77) = 2 for all rj. Our assumptions read 
as follows: 

(Al-2) ll > 0, A > -2/i; 
(A3) A + ji ^ 0; 
(A4) rV-(A + M)- 

Theorem 3.1 from [2] implies the well-known fact that the solution to the thermo-elastic system 
can be decomposed into three components, one related to the parabolic eigenvalue A + Iji and 
two related to the hyperbolic one. The hyperbolic components satisfy wave equations and the 
corresponding dispersive estimate with L^-LF decay rate 7j _1//2 . The parabolic component decays 
like t~ x . The decomposition is just the usual Helmholtz decomposition. 

Corollary. Isotropic media have the dispersive decay rate t~ 1 / 2 . 



4. Cubic media 

Cubic media are defined by n = t-i and o~\ = 02 = 0, i.e. 

K " \(X + jjL)rj im (r - ji)r, 2 2 +fij 1 ' 

This matrix is positive if and only if 

tr A{rj) = t + ji>0 (4.2) 

and 

det A( V ) = ((t — li) 2 — (A + ji) 2 )rj 2 r, 2 + ji(r - /i) + /i 2 > (4.3) 

for all 776 s 1 . Using < ?7 2 ?7| = cos 2 <^sin 2 <f> < j we see that this is equivalent to fir > and 
AjiT + (r - ji) 2 - (A + ji) 2 = t 2 + 2/ir - A 2 - 2fiX = (t - A)(r + A + 2/z) > 0. 
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Condition (A3) means that the eigenvalues of A(rf) are distinct. On the contrary, A(rj) has 
multiple eigenvalues for a direction r\ if and only if A{rj) is a multiple of the identity or 

(A + /iJjjiTja = 0, and (r - = (r - /i)^. (4.4) 

The first equation means (A + fx) =0 or 77x772 = 0. In the latter case the second equation is 
satisfied for (r — /z) = 0, in the first one for r\\ = r/J. These conditions describe degenerate 
directions and (A3) holds if and only if (A + /j)(t — /j) 7^ 0. 

In order to find all hyperbolic directions we have to solve 

(r - A — 2/z) sin 40 = 0. (4.5) 

The case r = A + 2[i gives isotropic media and we can skip it. In the remaining cases we 
require sin 40 = 0, i.e. = kit /A. This corresponds to eight hyperbolic directions and using 
the general statements we immediately know that the coupling functions vanish to first order in 
these directions. 

It remains to check condition (A4). We distinguish between the directions = 0, i.e. r) = (1,0) T , 
and = 7r/4, i.e. r\ = (V2/2,V2/2) T . In the first case A(rj) — diag(r, /j) and the hyperbolic 
eigenvalue (i.e. to the eigenvector rj 1 - = (0, 1)) is /1. Thus we require 7 2 ^ /i — r. For the second 
hyperbolic direction we obtain similarly the eigenvalues spec A(rj) = {/1 + (A + t) /2, (r — A) /2}. 
Again the second one corresponds to r/ 1 - and is therefore the hyperbolic eigenvalue. Thus we 
have to require "f 2 =/= (r - A)/2 - fj, — (A + r)/2 = — (/i + A). 

We collect our assumptions: 

(Al-2) [j,,t > 0, -r- 2/x < A < t; 

(A3) (A + M )(r-/i)^0; 

(A4) 7 2 ^{-(A + M ), /J -r}. 

Figure [l] depicts Fresnel curves and coupling functions for different cubic media. 
Corollary. Cubic media have the dispersive decay rate i -1 / 2 . 



5. Rhombic media 
Rhombic media are determined by the condition <j\ = 02 = 0, i.e. 

Ain) = ( {T }-/ ] f + ^ ( {X + ifT ) ■ 

W; V ( A + M)'7i ? 72 (t 2 - HM + fij v ; 

This matrix is positive if and only if (again using polar co-ordinates) 

tr A(rj) = (i + t 2 + (n - r 2 ) cos 2 > 0, (5.2) 
i.e. if fj, + min(ri, r 2 ) > 0, and 

det A(r?) = \it\ cos 4 + /jt 2 sin 4 — (A 2 + 2A/_t — tit 2 ) cos 2 sin 2 > 0, (5.3) 

i.e. if /jri,/jT2 > and 2/i^/riT 2 > A 2 + 2A/i — tit 2 . The last inequality can be simplified to 

(Vn^2 - A)(Vr[7a + A + 2/x) > 0. 

Assumption (A3) means that we find no solutions to 

(A + ji) sin 20 = and (ti — fi) cos 2 = (r 2 — /i) sin 2 0. (5-4) 

Thus the first condition gives that (A3) is violated ifA + /i = Oor0 = kn/2. In the latter case 
the second equation gives T\ = /1 or r 2 — fi (depending on which we chose) . If on the contrary 
A + /j = 0, the second equation reads as (ri + r 2 — 2/i) cos 2 = r 2 — fx, which has further solutions 
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Figure 1. Frcsncl curves for cubic media together with the corresponding cou- 
pling functions (as curves {(2 + aj(ri))r] \ r\ £ S 1 }). Depicted are A = 0, \x = 1, 
r = \ (upper left), A = 0, yU=l, r=| (upper right), A = —2, /j, = 1, r = \ 
(lower left) and A = —2, \x = 1, r = | (lower right). 
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if either t\, t 2 > \i or ti,t 2 < /U. Thus for (A3) we require that (n — ^)(t 2 — fi) ^ together 
with (A + n) ^ if (n - h)(t 2 - m) > 0. 

Hyperbolic directions are determined by 

2(n - r 2 )sin20+ (n + r 2 - 2A - 4^) sin 40 = 0. (5.5) 

This equation has the four zeros <f> — kir/2 and, dividing by sin 20 we see that all further 
zeros satisfy (n — r 2 ) + (n + t 2 — 2A — 4yu)cos20 = 0. This condition gives further zeros if 
t\ - t 2 < ti + r 2 - 2A - 4,u or t 2 - ti < n + r 2 - 2A - 4/z, i.e. if (A + 2/i - n)(A + 2jU - r 2 ) > 0. If 
this expression vanishes, but t\ ^ r 2 , we get triple zeros for the corresponding 0. If t\ = A + 2/i, 
T"2 7^ Ti this holds for = kir, otherwise for = kir + |. 




Figure 2. Fresnel curves and coupling functions for a rhombic medium with 
A = |, /z = 1, Ti=4 and r 2 = 8. The medium has eight hyperbolic directions. 




Figure 3. Fresnel curves and coupling functions for a rhombic medium with 
A = 1, fi = 1, Ti = | and r 2 = 8. The medium has four hyperbolic directions. 

To check condition (A4) we have to determine the corresponding hyperbolic eigenvalues. For 
= this gives 7 2 ^ fj, — Ti, for = J similarly 7 2 ^ /i — r 2 . If the additional direction 
= j arccos( Ti + 7 2 1 ~ 2 7 ^_ 4 ^ ) exists the difference of the eigenvalues is A + /i (as for cubic media!) 
and we require further j 2 ^ — (A + ^i). 

In detail: We know that the eigenvectors of A{rj) are given by rj = (cos 0, sin 0) and rj 1 - = 
(— sin 0, cos 0). This allows to determine the eigenvalues by just calculating the products A(ri)rj 
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and A^r/- 1 . We give only the first components, in the first case it is (7*1— /x) cos 3 <fi+[icos (f>+(\+ 
fi) sin 2 0cos0, which gives the (parabolic) eigenvalue (ri— /i) cos 2 0+(/i+A) sin 2 </>+/j = {t\ — 2/i— 
A) cos 2 0+A+2/i. For the second one we obtain (— (ri— /i) cos 2 ^>sin 0— /itsin</>+(A+^) cos 2 0sin</>, 
which by division through — sm4> gives the hyperbolic eigenvalue (r x — 2/i — A)cos 2 +/i. Thus 
the difference between both is — (A + /x). 



Figure 4. Fresnel curves and coupling functions for a rhombic medium with 
A = «, /x = 1, n = 1 and 72 = 8. The medium has four hyperbolic directions, 
in two the coupling function vanishes to third order. 

Except for the case with triple zeros, i.e. for (A + 2/i — ti)(A + 2/i — T2) = 0, we are done and 
the results from [5] imply the dispersive decay rate i" 1 / 2 . In the exceptional case we have to 
determine further the tangency index 7(77) of the corresponding sheet of the Fresnel curve. Since 
for rhombic media the Fresnel curves are mirror symmetric with respect to the co-ordinate axes 
and the exceptional directions coincide with one of these axes, we can only have 7/(77) = 2 or 
7,(77) = 4 for the corresponding sheet. Because the hyperbolic eigenvalue in the exceptional 
direction is /i and the parabolic one is r = A + 2/i this means that if A + /i < we are on the inner 
sheet and 72(77) =2. If A + /i > we are on the outer one and the higher order tangency occurs 
according to [2j Prop. 3.4] if d^y/ ' xi(4>) + \/ >c\ {4>) = for the corresponding 4>. We consider 
the case T\ = A + 2/x such that the coupling function 01(77) vanishes in <j) = to third order. To 
decide about the order of tangency it suffices to determine xi(<fi) up to 0(<fi 3 ). 

We use 77 = (cos</>, sin0) T = (1 - \<j) 2 ,(f>) T + 0(<p 3 ), such that 

^)^ (A+ (aU 2) + M (r!-t)^ + ,) + °^ (5 ' 6) 

nw = ( 1 rf 02 )+o(0 3 ), (5.7) 

xi(<j>) = /7 + /30 2 + O(0 3 ) (5.8) 
with an unknown constant (3 to be determined via A(0)r 1 ((/)) = >f 1 (0)r 1 (<^)+O(0 3 ). This implies 

v /i + (r 2 - A - \n)4> 



AO^HI , ,.A+0{4?) (5.9) 



and 

v /i +(/3- 2(l)(p 



M0)riW = [..^^1 , AA 2 ) + 0(4>") (5.10) 
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such that the unknown constants are determined by 

/ 3-I M = T2 -A-|/x. (5.11) 

This implies (3 = t 2 — A — 2/i. Thus xi(4>) = + (t2 — A — 2/t)0 2 + 0((j> 3 ) (even by symmetry of 
the Fresnel curve modulo 0(4> A )). We use that 

(0)) 2 



2 + v^)) 



= dix 1 (<f>)- 



•2xi(0) 



i=0 



5=0 



2xi (0) 

= ^Xi(O) + 2xi (0) = 2(r 2 — A — 2/i) + 2/i = 2(r 2 - A - fi). (5.12) 

This expression vanishes if and only if t 2 = A + /t. Hence, if r 2 =/= A + /i we get 71(77) = 2 in this 
direction, while otherwise 71(77) = 4 holds true. 




Figure 5. Fresnel curves and coupling functions for a rhombic medium with 
A=|,/x = l, ri = | and r 2 = §• The medium has four hyperbolic directions, 
in two the coupling function vanishes to third order and the Fresnel curve has a 
higher order tangency. 

Finally we collect the conditions we have to impose for rhombic media: 
(Al-2) M>n,r 2 >0, - A /rTT^-2/i< A< y^v^; 
(A3) (A + M)(Ti-/x)(r 2 -/i)^0; 

(A4) 7 2 £ {/i-Ti,/J-T 2 ,-(A + /z)} if (A + 2/t-n)(A + 2/i- r 2 ) > and 
7 2 g {/1 - n,/x - t 2 } if (A + 2/i - n)(A + 2/x - r 2 ) < 0. 

In Figures [2] and [3] correspond to the situation of simple zeros of the coupling function, while in 
Figures [4] and [5] one coupling function vanishes to order three. 

Corollary. Rhombic media have the dispersive decay rate i" 1 / 2 , except in the case t\ = A + 2/i 
and t 2 = A + /i (or vice versa), where the dispersive decay rate is t~ x l A . 



6. An exceptional type of media 

We want to consider media with t\ = t 2 = A + 2/i, <7± — and cj 2 = /t. Again we start to check 
conditions (Al) to (A4). The matrix ^(77) has the form 



A(n) 



(A + /J)r?f + ,u 
(A + n)i]ir/2 



ml 



(A + [i)r]\r\i 

(A + + /i + 2/7771772 



(6.1) 
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such that it is positive if and only if 

tr A{rj) = A + 3^ + /i sin 20 > 0, (6.2) 

i.e. \i < A + 3yU, and 

det A(ri) = + 8A + 4^ cos 20 - cos 40 + 8,(isin20 + 4(A + fi) sin 40) > 0. (6.3) 

8 

This is clearly positive if A > \i > 0. [One can relax this condition to A > — s ~^* n anc ! ^ > 0. 
Note, that -2 + §\/3 w -O.7.] 

Assumption (A3) is violated if and only if for some directions 

(A + ^)iliV2 + = and (A + /u)r?i = (\ + n)^ + 2fii] 1 i] 2 - (6.4) 

The first equation implies % = or (A + fi)r]i + ^n] 2 = 0. In the first case the second equation 
implies A+/i = 0, while in the second one the second equation implies 2(A+/j)?7 2 = A+/i+2/i77i?7 2 , 
which together with the first one gives = A + \i — 4(A + Thus either A + = or ?? 2 = \. 

This corresponds to rfe = I an d thus we need = A + fj, ± v^Mj i- c - A + (1 ± \/3)m = 0. The 
condition with + can be skipped (since it is excluded by (A2)), the same for A + fi. Thus for 
(A3) we require A ^ (\/3 — l)/i. 

Hyperbolic directions are determined by 

= 4^icos20 - 4yticos40, (6.5) 

i.e. cos20 = cos40 = 2cos 2 (20) — 1, which implies cos20 = j(l ± 3) £ {1,— |}. Thus the 
hyperbolic directions are = kir/3 and we got double zeros for = kir. We check condition (A4). 
For = (or — tt) the eigenvalues of the matrix A(rj) are /i (hyperbolic) and A + 2/x, thus we 
exclude 7 2 ^ — (X+fi) (which never happens since by (A2) A+^i > 0). For = 7r/3 (or = 47r/3) 
they are fj, — \/3/u/4 (hyperbolic) and A +2/x+ 3-^^/4, thus we exclude 7 2 7^ —fj, — X — \/3fi (which 
also never happens due to (A2)). For = 27r/3 (or = 5n/3) we get /i + \/3/i/4 (hyperbolic) 
and A + 2/i — 3^/3^/4, thus we have to exclude 7 2 7^ —A — /i + \/3[i. 

It remains to consider the exceptional directions = kir and to investigate the vanishing order 
of the coupling functions and the tangcncy index of the Fresnel curve. We have seen that the 
hyperbolic eigenvalue ji is always smaller than the other eigenvalue, thus we are on the outer 
sheet of the Fresnel curve. Similar to the treatment of rhombic media we determine >?i(0) modulo 
O(0 4 ). Because 77 = (1 - ±0 2 , - ±0 3 ) T + O(0 4 ) we get 

*M~*) + '(£) + *(£) +ow ' ) «»> 

x x (0) = fi + /?i0 + /? 2 2 + /3 3 3 + O(0 4 ) (6.8) 

with unknown constants ol\ , 0:2, 013, CX4 and Pi,p 2 ,p3- Using 

_ (-iuf> + (ai(A + 2/x) + ^)0 2 + (a 3 (A + 2 M ) + a 2 (A + /1) + |(A + ^))0 3 \ 4 
A W n(0)-^ Ai + 2 M + a2A1 2 + ( a4M + ai (A + M ) + 2 a2/ i-| M )0 3 J+^J. 

(6.9) 

„ aw al f -^+(Mai-/3i)0 2 + (^ 3 + /3iai-/3 2 )0 3 \ „, ,<k , fi 1fV . 

*i(0n(0 - ( + M + {m2 + m , + {mA + ftaa + ) + 0(4, ) (6.10) 
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Figure 6. Fresnel curves and coupling functions for the exceptional medium 
with A = and fi = 1. The medium has six hyperbolic directions, in two the 
coupling function vanishes to second order and the Fresnel curve has a point of 
inflection in these directions. 



we see that /3i = 2fi, a± — — 3^/(A + fi), fit = 0, (3^ = — 7/3 (while we would need an 
additional normalisation condition to determine 02, «3 and a^). Thus 

xx{4>) =[i + 2ii4,- (3/i + l)^ + O(0 4 ). (6.11) 



This allows to determine 71. Note first that 
Thus 71 (f?) > 2. On the contrary. 



= ^ l( 0)-^ffl + 2x l( Q) 
0=0 2x H°) 



2^ K^d^dl^M + v^iM) 



4=0 



02*l(O) + (fy*!(0))(l 



-2/i + 2// = 0. (6.12) 
35^(0) 3(^x1 (0)) 3 



2^(0) 



4(*i(Q)) 2 



24^ 3 
V 



-(18^ + 21) +2^+ = -10^- 21 < (6.13) 



and therefore 71 (rj) = 3 for this direction. 

Again we collect the conditions for this exceptional type of media 
(Al-2) A > fi > 0; 
(A3) A^(^-1)m; 
(A4) 7 V(x/3-l)M-A. 

In Figure [6]the Fresnel curves and the coupling functions are depicted for one exceptional medium 
of this kind. 

Corollary. This exceptional medium has the dispersive decay rate i^ 1 / 3 . 
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